#include<bits/stdc++.h>
#define INF 0x3f3f3f3f
using namespace std;
typedef long long LL;
typedef double db;
inline int read()
{
	int x=0,f=1;
	char ch=getchar();
	while(ch<'0'||ch>'9') { if(ch=='-')f=-1; ch=getchar(); }
	while(ch>='0'&&ch<='9') { x=(x<<1)+(x<<3)+ch-'0'; ch=getchar(); }
	return x*f;
}
const int MAXN=110,MOD=9999973;
int n,m;
LL f[MAXN][MAXN][MAXN],ans;
inline int C(int x)
{
	return x*(x-1)/2;
}
int main()
{
	n=read();m=read();
	f[0][0][0]=1;
	int cur=0;
	for(int i=1;i<=n;++i)
	{
		cur^=1;
		for(int j=0;j<=m;++j)
			for(int k=0;j+k<=m;++k)
			{
				f[cur][j][k]=f[cur^1][j][k];
				if(j)(f[cur][j][k]+=f[cur^1][j-1][k]*(m-j+1-k)%MOD)%=MOD;
				if(k)(f[cur][j][k]+=f[cur^1][j+1][k-1]*(j+1)%MOD)%=MOD;
				if(j>1)(f[cur][j][k]+=f[cur^1][j-2][k]*C(m-j+2-k)%MOD)%=MOD;
				if(k>1)(f[cur][j][k]+=f[cur^1][j+2][k-2]*C(j+2)%MOD)%=MOD;
				if(k)(f[cur][j][k]+=f[cur^1][j][k-1]*(m-j-k+1)%MOD*j%MOD)%=MOD;
			}
	}
	for(int i=0;i<=m;++i)
		for(int j=0;i+j<=m;++j)
			(ans+=f[cur][i][j])%=MOD;
	printf("%lld\n",ans);
	return 0;
}
